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'P' Abstract 

Oh! 

r-| ! It has recently been shown that a minimal SO(IO) model with a single 10 and a single 126 

Higgs field breaking B-L symmetry predicts large solar and atmospheric mixings in agreement 



with observations if it is assumed that the neutrino mass obeys the type II seesaw formula. No 
additional symmetries need to be assumed for this purpose. Understanding CP violation in the 
renormalizable version of the model however requires a significant non-CKM source. In this 
paper we show that if we extend the model by the inclusion of a heavy 120 dimensional Higgs 
field, then it can accommodate CKM CP violation while remaining predictive in the neutrino 
sector. Among the predictions are: (i) solar mixing angle in the observed range; (ii) 6*13 in the 
range of 0.1 to 0.26; (iii) the Dirac phase close to maximal for the central value of the solar 
mixing angle. 



1 Introduction 

The simplest grand unified model for understanding small neutrino masses appears to be the 
SO (10) model PP for the following reasons: (i) it automatically brings in the right handed 
neutrino, Nji, needed to implement the seesawj^I mechanism since it fits in with other stan- 
dard model fermions in the 16 dimensional spinor representation (ii) it contains the SU{4:)c 
symmetryjn] which relates the quark and lepton coupling parameters and in turn helps the 
predictivity of the model in the neutrino sector by reducing the number of parameters; (iii) it 
also contains the B-L symmetry [21 E] needed to keep the right handed neutrino masses below 
the Planck scale and provides a group theoretic explanation of why neutrinos are necessarily 
Majorana particles. 

While all these make the SO(IO) models appealing for neutrino mass studies, detailed quan- 
titative predictions generally involve too many parameters limiting the predictive power unless 
extra symmetries (e.g. family symmetries etc) are imposed on the theory. One exception to 
this is the class of models that uses only one 10 and one 126 Higgs multiplet to generate 
fermion massesjSj. The original set of papers on this model[Hl IE] used type I seesaw formula 
for neutrino mass is given by Adi, = — Mj^ M^^ (M^)^ where M^ is the Dirac mass of the 
neutrinos and Mr is the mass matrix of the right handed neutrinos. These predictions are now 
in contradiction with experiments. 

It was subsequently pointed out in ref. [7j that if one uses type II seesaw formula for the 
neutrino masses jS] instead, the model automatically predicts large atmospheric mixing angle 
due to the fact that bottom quark and tau lepton masses converge towards each other when 
extrapolated to the GUT scale. The question remained whether this works for three generations 
and can lead to a realistic model for neutrinos. It was shown in ref. [H] that the same b — r mass 
convergence not only leads to a large solar mixing angle, but also to a small and detectable 
value for Ue3 = sin ^13. A detailed numerical analysis was carried out that showed that the 
model is indeed in agreement with present neutrino data and in particular the prediction of a 
"large" value for 6*13 which makes this model testable at the current as well as at the proposed 
long base line neutrino experiments. 

In the three generation neutrino discussion in ref. jH] , the Yukawa couplings of fermions were 
assumed to be real and all CP violating effects were assumed to originate from the supersymme- 
try breaking sector. It is however interesting to check if one can accommodate the CKM phase 
in the model by introducing phases in the couplings. A detailed investigation of the minimal 
model where CP violation is introduced through complex Yukawa couplings (as in the standard 
model) showed PU] that compatibility with neutrino data requires the CKM phase to be outside 
the first quadrant whereas the standard model CKM phase is in the first quadrant fT]. This 
would seem to imply that in order to understand observed CP violation in this model, one 
must invoke a significant non-CKM source for CP violation (as in the model with real Yukawa 
couplings) e.g. CP violation from the supersymmetry breaking sector. This could very well be 
true. However, since all observed CP violating phenomena seem to be explainable by the CKM 
model, it is important to see whether one can explain both CKM CP violation and neutrino 
mixings by a minimal modification of this SO(IO) model. There are also other issues such as 
SUSY CP problem that one needs to address in the context of supersymmetry and it would be 
interesting to see how these can be addressed in this model. 

In this paper, we propose a very minimal way to incorporate CP violation into the model. 



which not only leads to a predictive model in the neutrino sector but also seems to have wider 
implication beyond just explaining CKM CP violation. For instance, the model presents a 
solution to the SUSY CP problem. 

In order to attain our goal, we include a heavy 120 field with an extra Z2 symmetry which 
we will call "parity" symmetry imposed on the model^. At energy scales below the mass of 
the 120 field, the effect of this field is to appear as a higher dimensional contribution to the 
Yukawa couplings. This effective theory has the following properties. Despite the fact there 
are now three extra parameters in the model, the theory still remains predictive in the neutrino 
sector. Secondly, the mass matrices for quarks and leptons are hermitian, which therefore has 
the potential to solve other problems of supersymmetric models such as the SUSY CP problem. 
In this paper we focus only on the neutrino sector. 

The main results of this paper are as follows: (i) using type II seesaw formula we are able 
to accommodate the CKM CP phase while keeping the model predictive in the neutrino sector; 
for example, we predict the solar mixing angle in the right range and Vez > 0.1; (ii) the model 
has the potential to solve the SUSY CP problem and (iii) it predicts the Dirac phase of PMNS 
matrix to be near maximal for the central value of the solar mixing angle tan^^0 ~ 0.4. 

The paper is organized as follows: in sec. 2, we introduce the model with the inclusion of 
the 120 Higgs field and write down the fermion mass formulae in the general case; in sec. 3, 
we impose parity symmetry on the model making it predictive in the neutrino sector; in sec. 
4, we discuss the predictions for neutrino mixings and Dirac CP phase for neutrinos; in sec. 5, 
we present our conclusions and discuss the outlook for the model. 

2 SO(IO) Model and CP violation 

We start by writing down the Yukawa interactions of our model, which are responsible for the 
discussion of neutrino masses and mixings. The Yukawa superpotential involves the couplings 
of the 16-dimensional matter spinor ^j with 10 (i/), 126 (A), and 120(A) dimensional Higgs 

fields: 

1 1 _ 1 

Wy = -hijtpitfjjH + -f.jtpitpjA + -h[jiji^ljjA. (1) 

The Yukawa couplings, h and /, are symmetric matrices, whereas h' is an anti-symmetric matrix 
due to SO(IO) symmetry. They are all complex matrices in general. 

Once the SO(IO) symmetry breaks down to the standard model symmetry, we have four 
pairs of Higgs doublets arising from the H, A, and A Higgs fields. There may also be other Higgs 
doublets e.g. in 210 multiplet. Under the G422 = SU{4:)c x SU{2)l x SU{2)r decomposition we 
have the following representations that contain the Higgs doublets of up and down type: one 
pair arises from H D (1, 2, 2), one pair comes from A D (15, 2, 2), and two pairs come form 
A D (1, 2, 2) + (15, 2, 2). We assume that one pair of their linear combinations, Hu and Hd-, 
remains massless (mass is ~ Oivwk)) and become the MSSM Higgs doublets. As for other pairs, 
they all have GUT scale masses. Using the light Higgs doublets, the MSSM Yukawa couplings 
below the GUT scale and the right-handed Majorana neutrino mass terms can be written as 

Wy D Yi;Q,U'rHu + Y.JQ.D'^Hd + Y^LiE^H^ + Y.^^UN'rH^ + ^/.,L,L, A^ + ^/..iV^^iV^A^, (2) 



^Vcry different extensions of the model that use 120 but no symmetry have been discussed in ref.^S]- 



where Q, If^, V^, L, E'^, N"^ are the quark and lepton superfields which are all unified to the 16 
spinor tp field. A^ is an SU{2)l triplet Higgs field and /S.j^ is a neutral component of SU{2)ji 
triplet, both part of the 126 field. Even though both of them have GUT scale masses we have 
included them with the MSSM superpotential because their vevs lead to light neutrino masses 
via the seesaw mechanism. 



The gauge coupling unification requires that the 126 Higgs field acquires VEV at or close 
to the GUT scale. We also need to introduce a 126 Higgs field to satisfy the D-fiat condition to 
maintain supersymmetry down to the weak scale. Though the 126 Higgs field does not couple 
to the fermions, the pair of Higgs doublets in the 126 mix with the doublets arising from the 
other Higgs multiplets, since 126 couples to the other Higgs multiplets with non-zero coupling. 
These five pairs of Higgs doublets, H]^^, Au^, A^^^, A^ ^ and A" ^ are mixed and the light pair 
of Higgs doublet can be written as 

(i7„,---) = (i7^°,A„,A„,A:,Af ,■■■) Uh, (3) 

(i/,,---) = (iJf,A,,A,,A^,Af ,■■■) Vh. (4) 

where Uh and Vh are unitary matrices. The superscripts s and adj stand for SU{'i)c singlet and 
adjoint pieces. We have temporarily ignored the doublets that may arise from other multiplets 
in the theory such as 210. It is important to stress that in order to obtain one pair of MSSM 
Higgs doublets from five pairs at GUT scale, one needs to do a fine tuning of parameters. We 
have enough parameters in the Higgs superpotential that this is possible to achieve. We have 
also checked that we do not have any light color triplet fields. 

The results given below remain unchanged in their presence. The Dirac mass matrices of 
quark and leptons are 

Mu = Mio + Mi26 + Mi2o, (5) 

Mrf = riMio + r2Mi26 + rsMi2o, (6) 

Me = riMio - 3r2Mi26 + Ar^Muo, (7) 

M^ = Mio - 3Mi26 + AM120, (8) 

where the three mass matrices in the expression are given by Mio = h*Vu{UH)u, M126 = 
Cif*Vu{UH)i2, and M120 = h'*Vu{{UH)ii + c-2{UH)ib) where Vu is a vacuum expectation value of 
MSSM Higgs doublet H^-, and Cj are Clebsch-Gordon (CG) coefficients. The coefficients r^ and 
A are written as 

r-i = {VH)ii/{UH)iicotf3, (9) 

r2 = {VH)lz/{UH)l2COi(3, (10) 

ra = {{Vh)ia + C2{VhU)/{{Uh)ia + C2{UHU)coil3, (11) 

r4 = {{Vh)ia - ^C2{Vh)i^)/{{Uh)ia - 3c2(t//f)i5) cot/5, (12) 

A = {{UhU - ^C2{Uh)i,)/{{Uh)ia + C2{UhU), (13) 

where cot/3 is a ratio of vacuum expectation values of doublet Higgs fields, cot/3 = Vd/vu- The 
Majorana mass matrices of left- and right-handed neutrinos prior to seesaw diagonalization are 
given by 

ML = rvL. MR = rvn, (14) 



where vl and vr are vacuum expectation values of Al and A/j, respectively. As already men- 
tioned, since vr is expected to be close to the GUT scale, this implies that vl is ~'y^cak/('7^Gf/T) 
<C fweak; wherc T] IS B. coupllug coustaut in the Higgs potential. The Majorana mass matrix of 
the heavy right handed neutrino is proportional to Mi26- The light neutrino mass matrix is 
given by the mixed type II seesaw formula, 

A^«sht_jVfi-M,^M^^(M,^)^. (15) 

As discussed in earlier papersjUinj, there are regions of the parameter space in the theory where 
the first term will dominate; we will call this the pure type II seesaw case. If on the other hand, 
we consider the parameter space where the second term is dominant we will call this type I 
seesaw. The bulk of our results will be for the pure type II case. 

3 Parity invariance and a predictive model for neutrinos 

In order to see if the model is predictive for neutrinos, let us count the number of parameters in 
the theory. In the basis, where M126 is real and diagonal, there are 3 real parameters in M126, 
6 complex parameters in Mio and 3 complex parameters in M120. We also have 5 complex 
parameters in the Eqs.ipiT^ as well as the vevs of the Altogether, there are 31 real parameters 
in the fermion sector, and, therefore, we do not have any prediction for the neutrino mixings. 

In order to be predictive in the leptonic sector of the model without imposing any flavor 
symmetry, we require the theory to be invariant under a parity symmetry. As we will see, it 
makes the Dirac mass matrices Eqs.(jnilHl) hermitian and leaves a total of 17 real parameters in 
the fermion sector making the model predictive. If we further require that the 120 Higgs field 
has a mass much higher than the GUT scale, its only manifestation is as an effective dimension 
four term in the superpotential. The reduces the number of parameters to 15 increasing the 
predictive power of the model. We explore both the cases with 17 and 15 parameters in a 
subsequent section. 

We now define the parity transformation in the G422 basis. We write SU(4) indices by /i, u, 
SU(2)i indices by a, [3 and SU(2)i:j indices by d, (3. The SO(IO) spinor ifj and x ^-^e decomposed 
as 

ij = ijf,a + ij^, X = X,ia + XA- (16) 



Bi-doublet Higgs fields in the 10, 126 and 120 are written as Haa, A^^aa, Aaa, and A^'^aa- 
Then the Yukawa interactions are written in the following (up to overall factors) 

h iJxH = h{^,^x'a + raX,c.)ir^ + ■■■ (17) 

/ V^xA = /(^M-X ■ + CXMa)A/"" + ■ ■ ■ (18) 

h' i,xA = h\,lj^^xl - CX/.a)^"° + C2h'{,lj,^xl - CXm«)^/"" + ■ ■ ■ (19) 

where C2 is a CG coefficient. The Lagrangian is written as 



c= I d'ew+ I d'ew (20) 

and 

c= fd'e h{4j,^x^ + r^x,a)H^'' + / d''e h*{{4j,^r{x^r + m*{x,ar){Hn* + ■■■ m 



We consider the symmetry under the following parity transformation, 

^^M« ^ ivn*, d'e ^ (fe. (22) 

Of course, x is also transformed in same manner. 

In the Higgs sector, the transformations of the (1, 2, 2) and (15, 2, 2) sub-multiplets under 
6*422 are: 

A consequence of the parity symmetry (J2S1); is that the coupling matrices h and / real and 
symmetric and h' antisymmetric and imaginary; the parameters ri{i = 1,2,3,4) and A in the 
Eqs.ll MT!^ are real. This considerably reduces the number of parameters in the theory and 
further makes the mass matrices for all charged fermions hermitian. 

Let us clarify our motivation for introducing the 120 Higgs field. When the 120 Higgs field 
is absent, the fermion mass matrices are complex symmetric matrices in the absence of the 
parity symmetry and we have the following relation in the pure type II case 

M,(xMd- nMu = UiVDdV^ - nDuW 



U 



\ Vuhtrib Vcbrrib nib - rifUt J 



U' , (24) 



where rric and rriu contributions in (2,2) and (1,1) elements are omitted, and (pd and (ps are 
complex phases in the diagonal matrix, Dd- If Me is close to a diagonal matrix in the basis 
where M^ is diagonal, the maximal atmospheric mixing can be easily obtained when the (3,3) 
element is suppressed such that Im^e*''^'' — (m(, — rimt)\ ^ 2Vcb'mb. This suppression of (3,3) 
element is related with other observed facts such as bottom and tau mass convergence at GUT 
scale and Aml^i/Am^ ^ 0{m'j./ml). Assuming that the atmospheric mixing is maximal, we 
obtain the neutrino mass matrix, Eq. ()24j) . as 



/ mde'^^ + Vlm^e'^^ -(Kb - KsKfe)mb/y2 (Kb - KsKb)mfe/v^ \ 



UU, 
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UlsU^, 



-(Kfe - VusVcb)mb/V^ evfiQ 

V (Kb - VusVcb)mb/^ mo j 

(25) 
where mo and emg are eigenvalues of (2-3) block, and e^ ~ Aml^J Ain\ and rriQ ~ 2'ms. Thus, 
the solar mixing and 13 mixing are proportional to |Kfe/Kb — Ks| and tan26'so« ~ tan26'i3/e. 
Therefore, those mixing angles also depend on the KM phase, 5km- Since Kb = |Kb|e~*'^'^'^, 
the KM phase in the first quadrant gives a smaller value for solar mixing angle rather than 
in the second quadrant. In order to obtain the proper value of solar mixing angle, we have to 
choose a smaller value of e. However, in the model without 120, the e parameter, which is a 
function of the strange quark mass, is constrained due to the fitting of three charged-lepton 
masses (especially electron mass), and we do not have proper fitting of the solar mixing angle 
data in the case where the KM phase is in the first quadrant. We can verify the situation in a 
precise analysis in the pure type II case jTU]. In the type I case, things are more complicated, 
but it has been shown that it is not possible to fit the neutrino oscillation data in this model 



with the above mininial Higgs choice ^S]- The mass squared ratio is constrained due to the 
charged-lepton mass fitting and it cannot be small enough when the KM phase is in the first 
quadrant. The mass squared ratio is a free parameter in the model with 120 since the additional 
parameter A in the sum rules Eqs.(jS]|Hl) can fit the electron mass, and therefore we can explain a 
smaller KM phase. Thus, we employ the 120 Higgs field to explain the large solar mixing angle 
along with the KM phase in the first quadrant. Interestingly, even though we have introduced 
a new Higgs field, the number of parameters is less than the minimal Higgs choice with most 
general CP phases due to the constraint of parity symmetry. 

We further note that if the 120 field is heavier than the GUT scale, its effect on the physics 
at the GUT scale comes from a higher dimensional operator of the form ipTTTijjH^/M, where 
$ is a 210 Higgs field, so that we get the relation ri = r^ = r^. In other words, the choice 
ri = rs = r^ is not an adhoc choice but can be guaranteed in a natural manner. 

We now note that in the presence of the parity symmetry, since all the mass matrices 
are hermitian and also the /i-term and the gluino masses are real, the most dangerous graphs 
contributing to large electric dipole moment of the neutron are absent [TH]. Therefore this model 
has the potential to solve the so-called SUSY CP problem. 

We also wish to recognize that the Z2 CP symmetry we impose is broken at the GUT scale 
by the VEVs of 126, 210 and 45 Higgs fields which break SO(IO) symmetry. The light MSSM 
doublet Higgs fields are no more CP eigenstates, and thus there is no cosmological domain wall 
problem at weak scale in the model. 

4 Near Bi-Maximal Solution for Neutrino Oscillation 

As already noted in section 3, under the parity symmetry, Miq and M126 are real symmetric 
matrices, M120 is a pure imaginary anti-symmetric matrix, and the coefficients r^ and A are real 
parameters in the Eqs.(j3]|Hl). ^^ can therefore rewrite the charged-lepton and Dirac neutrino 
mass matrices in terms of the other mass matrices as follows: 

Me = Cu ReMu + Cd ReMd + lAr^ ImM„, (26) 

Q I 

M^ = ReM„ + Re(Md - nM^) + iA ImM„, (27) 

n 

where c„/(l — q) = ri, —Cu/{3 + q) = r2. The up- and down-type quark mass matrices are 
hermitian, and are written as 

M„ = f/A.f/^ Md = UVDdV^U\ (28) 

where Du = diag(±m„, ±mc,mt), Dd = diag(±md, im^, m^) and V is the CKM matrix, and 
f/ is a unitary matrix. We note that nit and rrib component in the Du and Dd can be made to 
be positive without loss of generality. Because of the parity symmetry, Md — r^M^ must be a 
real symmetric matrix. We fix the flavor basis as 

Md - rsM, = U{VDdV^ - r^D^W^ = diag. (29) 

The unitary matrix U is determined by r^, up to phase matrix P, 

U = PU, P = diag(e^'^^ , e'^^ , 1) (30) 



The parameters are now 6 quark masses (with signatures), 3 mixing and 1 KM phase in the 
CKM matrix, and the coefficients q, c„, 0i, 02? ^3, ^a, ^- There are 17 parameters in all. For 
example, the three charged-lepton masses can be used to determine q, c„, At/^. The remaining 
4 parameters give the neutrino oscillation parameters. 

In the case where the 120 Higgs field is heavier than the GUT scale, it manifests itself as an 
effective dimension four operator of the form ifjifjll^^ / M . As a result, the Higgs doublets in 120 
are decoupled. This leads to a reduction in the number of mixing parameters. This translates 
into the the relation ri= r^ = r^ since the vev ratio that defines t'i_3^4 is the same. In this case, 
there remain only two parameters describing the neutrino sector. They can be determined by 
two of the parameters from the neutrino oscillation data and the remaining neutrino parameters 
can then be predicted. 

Interestingly, if we have the relation ri = r^, the matrix Ma — r^Mu is proportional to the 
light neutrino mass matrix M^ and the diagonalizing matrix U become close to MNSP matrix 
in the pure type II case. Therefore, if ri = rs ~ mb/mt and the (3,3) element of Md — riM^ 
is suppressed, we have a large atmospheric mixing. The mass squared ratio is of the order 
of 10~^, which is the right order seen in the experiment, only if ri ~ mb/rrit (otherwise, the 
mass squared ratio become 0{m1/mf) or 0{m1/ml)). Furthermore, since the (3,3) element of 
Mi26 is suppressed for that choice of ri, the bottom-tau mass unification is satisfied and it is 
consistent with the renormalization group fiow for the case of tan/3 ~ 50. 

Now let us study the prediction of the model in the case where ri = r^ = r4^. In this case, 
we have 15 parameters in the model. After fixing the quark masses (with signatures) and the 
CKM parameters, we are left with 5 parameters, q. A, 0i, 02, and ri. Since mass squared 
ratio is a function of ri, we can fix the parameter ri by the experimental value of Aml^i/Am^. 
The three charged-lepton masses can be used to fix q. A, and 02- As a result, the neutrino 
oscillation parameters, 9a, Osoi, \Ue3\, and one CP phase Smnsp are predicted by only one phase 
parameter 0i. Interestingly, the atmospheric mixing 6a doesn't depend on the phase 0i very 
much, and the 6a is really predicted when we fix the quark masses and mass squared ratio of 
light neutrino. 

It should also be noted that the other arbitrariness in the model is due to the choice of the 
signs of different fermion masses, since the sign of a fermion mass is unobservable. We find 
that only for the two choices of the signs given below, we obtain acceptable solutions: 



(a) 


Du 


= diag(±, - 


-,+), 


Dd = diag(-,+,+), 


De 


= diag(±, -,+), 


(31) 


(b) 


Du 


= diag(±, - 


-,+), 


Dd = diag(+,+,+), 


De 


= diag(±,+,+). 


(32) 



The solutions we present correspond only for these two choices of signatures. 

In Fig.l, we show the prediction of the atmospheric mixing sin^ 2^a as a function of the 
mass squared ratio. The lines (a) and (b) in the figure are for the set of quark and lepton 
mass signatures given above which give acceptable solutions for neutrinos. We can obtain 
large atmospheric mixing angles with the proper choice of quark masses in their allowed range, 
mixings and KM phase ^Hl with the choice of the set of signatures, especially for case (a). The 
most important input parameter for obtaining a large atmospheric mixing is the strange quark 
mass. We show the strange mass dependence of 6a in Fig. 2 in the case where the mass squared 
ratio R = Amg/Am^ is 0.03. The strange mass in the figure is the running mass at 1 GeV. In 
order to obtain the experimental constraint sin^ 26a > 0.9, we need the parameter region where 



the strange mass has a larger value. The bottom quark mass dependence is not negligible, and 
a larger value of bottom mass is preferred to obtain maximal atmospheric mixing. 

After fixing all quark and lepton data and also the mass squared ratio, i?, we can fit the solar 
mixing angle by choosing the free phase parameter 0i. Then, \Ue3\ and Smnsp are predicted. 
In Fig. 3, we show the correlation between solar mixing tan^ 6soi and {Uesl by varying the phase 
parameter 0i. We give two lines for the cases (a) and (b), and the two lines correspond to 
different mass squared ratios, R = 0.02 and R = 0.07. From the figure, we can see that 
AmQ/Am^ = 0.07 is not favored in the 3a range of the experimental data of solar neutrino 
and Ue3 for the case (a). In Fig.4, we present the prediction of the sin Smnsp in the case where 
R = ArriQ/Am^ = 0.02. The CP phase can be of any value in the range of the experimental 
data of solar neutrino. If we restrict the mass signatures to the (a) case, sin 6mnsp could be 
predicted to be ±0.9. 

The most interesting feature of this model is the prediction of jf/esl- In Fig. 5, we show 
the bottom quark mass dependence of the jf/esl prediction. The bottom mass is defined as a 
running mass at rrib. Since a large SUSY correction to the bottom quark mass can be induced 
in the large tan/5 case, the running bottom mass can be large, and the larger bottom mass 
gives smaller value of jf/esl- In a similar way, larger tan/3 predicts smaller If/gsl since the larger 
tan/3 gives larger bottom quark mass at GUT scale. In Fig.6, we show the plot of {Uesl- Each 
point is dotted for different quark mass and mixings which are randomly generated in the 
experimentally allowed region. We can see that the model predicts a lower limit for l^esl of 
about 0.1. 

In Fig. 7, we can see the KM phase dependence of the prediction of the model. The lines 
in the figure are drawn by changing 0i for different values of KM phases. The mass squared 
ration is set to be 0.02. For a smaller KM phase, the lines shift to the smaller solar mixing. In 
this model, the experimentally allowed solar mixing can be obtained in the first quadrant KM 
phase, contrary to the minimal model without the 120 Higgs field. We can see that a smaller 
value of {Uesl can be allowed for the larger values of KM phase. 

Since all the parameters of the model are now determined, it can be used to make other 
predictions. As an example, we have calculated the magnitude of lepton flavor violating process 
/i — > e + 7 in the model (see Fig|H)). Note that the predictions are in the range currently being 
probed by experiments|21j. We use the mSUGRA universal boundary conditions at the GUT 
scale i.e. rrio (universal scalar mass), mi/2 (universal gaugino mass),/lo (universal trilinear mass 
). The other two parameters are sign of /i and tan/3. The dots in the plots are produced for 
various model points generated by fitting the fermion masses and mixing angles. We can see 
that the BR of /i — > 67 can be large for smaller values of mi/2- The lightest neutralino is the 
dark matter candidate in this model and we satisfy the 2a range of the recent relic density 
constraint I^cdm = 0.1126to;oo9|22] in the parameter space. When we satisfy the relic density 
constraint, the mo gets determined. We choose Aq = and n > 0. The right handed masses 
have hierarchies and therefore get decoupled at different scales. The flavor-violating pieces 
present in Y^, and / induces flavor violations into the charged lepton couplings and into the soft 
SUSY breaking masses e.g. m^ terms etc. Also, an additional symmetry between the GUT 
and the vr scale (type I) helps to induce flavor violation^^l- The electric dipole moment of 
electron is smaller than the experimental reach in the range of parameter space showed in the 
figure. 

We comment that in the case of type I seesaw, the large mixing solution is a sharp resonance 

9 



solution and the solution is not stable to predict the mixings contrary to the type II seesaw. 

5 Discussion 

In this paper we consider the prediction of neutrino masses and mixings for an SO(IO) model 
where fermion masses receive contribution from the presence of three Higgs multiplets 10, 120 
and 126. We impose a parity symmetry on the model, so that it has very few parameters which 
enables prediction of two mixing angles and all the CP phases in the neutrino mixing matrix. 
The advantage of this model over the most minimal SO(IO) model is that now the CKM phase 
is in the first quadrant as required by the standard model analysis of all observed hadronic CP 
violation. We also wish to emphasize that this one of the few models in the literature that can 
predict leptonic CP phases. 

As far as experimental tests of this model are concerned, the parameter f/gs = sin ^13 is 
predicted to be large like the most minimal SO(IO) modeljTl E] but is somewhat smaller i.e. 
Ue3 > 0.1. This can be tested in the next round of planned long baseline experiments. We 
also predict that the Dirac phase for neutrinos can be maximal. Furthermore the model has 
the potential to solve SUSY CP problem due to the fact that all fermion masses are hermitian. 
The model also predicts observable amount of LFV in muon decay. 
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Figure 1: The atmospheric mixing angle is plotted as a function of mass squared ratio. Predic- 
tions for different set of mass signatures (a) and (b) in Eqs. ()31ll3^ are given. 
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Figure 2: The atmospheric mixing angle is plotted as a function of strange quark mass. The 
strange quark mass is given as a running mass at 1 GeV. The mass squared ratio is 0.03 in each 
set of mass signature, (a) and (b). 



13 



0.30 



0.25 



0.20 



00 

-^ 0.15 



0.10 



0.05 



0.00 



~ 






- 


y 




- 


/ 




- 


/ 


^ — - — ■ ' 


■ 


/ / 


^ 


. 


/ 


^' 


/ 
/ 
. / 






- / 


I f f 






' J ■' 


/-l^ \rr? f \rr? — (\ HT 




" 




(a)^m'^/zlm^ = 0.07 


■ 


/ 


(b)^^,/dm^ = 0.02 


■ 


/ 


(b) zlm'^/zlm^ = 0.07 


■ 











0.0 



0.2 



0.4 



0.6 



1.0 



tan^0sol 



Figure 3: The relation between solar mixing angle and jf/esl is plotted. Each line is plotted 
by varing a free phase parameter 0i. The experimentally allowed region in 3a of recent data 
fitting is 0.3 < tan^ ^.^^ < 0.6 and \Ue3\ < 0.26 [20]. 
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Figure 4: The prediction of MNSP phase is plotted as a function of the solar mixing angle. 
These lines (a) and (b) are plotted in the case Am1^i/Am\ = 0.02 for different mass signature. 
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Figure 5: The relation between solar mixing angle and {Uesl is plotted as dots for different 
bottom quark masses. The bottom quark mass is given as a running mass at mi,. This plot is 
given in the case of mass signature (a) and mass squared ratio is 0.02. 
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Figure 6: The prediction of {Uesl ^^^ solar mixing angle is plotted as dots for randomly gen- 
erated quark masses (with signature) and mixings in the experimantally allowed region. The 
lower bound of It/esl exists in this model. 
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Figure 7: The relation between solar mixing angle and {Uesl is plotted for various KM phases. 
Each plot is given in the case of mass signature (a) and mass squared ratio is 0.02. 
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Figure 8: The BR[/i -^ 67] is plotted as a function of mi/2 for tan/3 = 10 in type II. 
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